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Abstract. We analyze a system of two qubits embedded in two different environments. The qubits are 
coupled to each other and driven on-resonance by two external classical sources. In the secular limit, we 
obtain exact analytical results for the evolution of the system for several classes of two-qubit mixed initial 
states. For Werner states we show that the decay of entanglement does not depend on coupling. For other 
initial states with "X"-type density matrices we find that the sudden death time displays a rich dependence 
on the coupling energy and state parameters due to the existence of processes of delayed sudden birth of 
entanglement. 
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In the theory of open quantum systems, when a two-level 
system (a qubit) is placed in a dissipative environment 
(such as a bath of harmonic oscillators), a number of deco- 
herence effects are known to occur and the mathematical 
formalism describing them has been worked out in the last 
few decades [T]. A very interesting situation occurs when 
one considers two qubits and their entangling properties: 
if the qubits are coupled to two different (uncorrelated) 
baths, qualitatively new effects, known as sudden death 
of entanglement phenomena, have been found to occur [2i 
[31[3]. These effects, already confirmed experimentally [S], 
include the non-exponential decay of concurrence and 
the nonadditivity of decay rates [3]. It has been shown 
that these sudden death effects becomes more pronounced 
at finite temperature [6] and in the presence of external 
driving [7]. Such systems have attracted considerable in- 
terest amongst various research communities, due to the 
fact that the entangling effects predicted are fundamental 
and generic - in the sense that they do not depend on a 
specific experimental realization. 

Much of the nonituitive character of the sudden death 
of entanglement comes from the fact that the system is bi- 
partite and, although it is subject only to local interactions 
(acting separately on each of the two subsystems), one 
property (entanglement) which characterizes the whole 
system, has a dynamics which cannot be understood qual- 
itatively as just the addition of local effects. There are 
also practical consequences to the sudden death effect: for 
the research effort into building a quantum computer, it 
implies that it could be more difhcult to maintain entan- 
glement when the qubits interact with different reservoirs 
and are locally driven by single-qubit quantum gates. In 
this context, also qubit-qubit interactions must be consid- 



ered: in quantum networks such as those realized nowa- 
days with superconducting quantum circuits these interac- 
tions can appear while performing two-qubit gates or sim- 
ply by spurious cross-coupling. This qubit-qubit coupling 
is established either directly, through mutual inductances 
and capacitances [5], or as second-order effects, due to the 
creation of virtual excitations in the common environment 
or in the quantum buses needed in a typical architecture of 
a quantum processor [S] . Finally, in the presence of driving 
and of two different reservoirs, the qubit-qubit interaction 
leads to another entanglement effect called steady-state 
generation of entanglement, a phenomenon which has at- 
tracted considerable attention [TD]. In contrast, transient 
entanglement effects that occur as shorter timescales in 
such systems have not yet been studied systematically. 
Much like in the case of steady-state entanglement, we 
will see that the interesting effects occur only when driv- 
ing, coupling, and decoherence are present at the same 
time. However, we should mention that the steady-state 
generation of entanglement appears in a completely dif- 
ferent region of the parameter space, in which the secular 
approximation is not valid: for the parameter values con- 
sidered here, the steady state is not entangled. 

In this paper, we investigate a system consisting of two 
qubits, dissipating in two different reservoirs, driven at 
resonance, and coupled by a dipole-dipole interaction. We 
study the entanglement dynamics of initial states of "X"- 
type, that is, 4 x 4 density matrices in which the nonzero 
elements have the shape of the letter "X". As we will see, 
under unitary evolution and in the secular approximation 
these states are mapped also into X-type states and find- 
ing analytical solutions is possible. For the states investi- 
gated, the dynamics occurs at timescales below the inverse 
of the longitudinal decoherence rate - the states become 
separable after a time of sudden death tESD- Bipartite 
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decoherence is truly a requirement: for a single reservoir, 
the decay would be exponential, which would make any 
diagram of tESD featureless. As it turns out, the three 
processes present in the problem - driving, coupling, and 
decoherence - combine into a non-trivial time evolution 
in the interval t G [OjIesd] resulting in an unexpectedly 
rich dependence of tESD on the initial state and on the 
coupling strength. 

We consider two qubits of energy separation (Larmor 
frequency) ajj" , and driven by two classical radiation fields 
of frequencies ujj which are coupled to the (T^-Pauli matri- 
ces of the qubits with dipole-field strengths (Rabi frequen- 
cies) ilj . The qubits are also dipole-dipole coupled to each 
other, a generic interaction which appears as such in many 
experimentally-relevant contexts [8,9J. In the Schrodinger 
picture we write for this part of the Hamiltonian 

H,oup^io-''<7f<7'^ + ioyy<7fai (1) 

where and cj^^ denotes the X-X and Y-Y coupling 
strengths, respectively. The total Hamiltonian of the sys- 
tem is therefore 



Here 6i = 



Uj is the detuning of the qubit transi- 



Hsys 



E 



i=i,2 



He 



(2) 



For the environment, we consider a longitudinal (am- 
plitude damping) reservoir at zero temperature. Under the 
usual assumptions of Markovian evolution with respect to 
the environment degrees of freedom [T], and at zero tem- 
perature, the dynamical evolution of the qubits' density 
operator (in the Schrodinger picture) ps is governed by 
the Born-Markov master equation 



PS = -i[Hsys,Ps] + C[ps], 



(3) 



where the resulting Liouvillean superoperator C is the sum 
of the longitudinal dampings of the two qubits with energy 
relaxation rates Fj , respectively, 
^ p. 

^[PS] = E ^ {^''JPS^'t - '^t'^JPS - PSCT+CT^) (4) 

We note here that we have implicitly used the typ- 
ical assumption (largely satisfied for many experimental 
systems) that the energy relaxation rates Fj have a mild, 
negligible variation with frequency over a bandwidth of 
the order of f2j around the Larmor frequencies of the 
qubits (or, leading to the same conclusion, one assumes 
that driving is weak with respect to the Larmor frequen- 
cies). We then move to a rotating frame using the trans- 
formation R = eicp[i(iOial -\- u;2cr|)i/2] and the identities 
R)a^R = eiq){^iu!jt) with j = 1, 2, and also employ- 
ing the rotating-wave approximation to eliminate the rel- 
atively faster counter-rotating terms, 

Prf = -i[Hrf, Prf] + C[Prf], (5) 

with Prf — RpsS^ , and 

.dR 



Hrf — RHsysR'' 



(6) 



E fl'^l + V^.' ) + T ^2 + <^M) ■ (7) 



tion frequency from the driving frequency and ujc is the 
effective coupling frequency (in the rotating wave approx- 
imation) LjJc = U!^^ + OJ^y ■ 

Eq. ([7]) can be solved numerically and we will present 
such results below. As it turns out, the data from nu- 
merical simulations will display qualitative differences for 
different initial states, and just the numerical result would 
be unsatisfactory. It is in principle possible to solve ana- 
lytically Eq. ([7]): unfortunately the structure of even rel- 
atively general forms of the density matrices (such as the 
"X"-form, see below) is not left invariant under evolution 
with Eq. ([7l), and the results are very complicated. To 
make progress, further simplifications are necessary. Be- 
low wc will use the secular approximation, which is valid 
for most of the interesting experimental situations when 
typically the Rabi frequency is much larger than the en- 
ergy relaxation rate, and we show that relatively easy to 
handle analytical results can be derived. 



2 Secular approximation for driven, coupled 
qubits 

In the rest of this paper we will analyze the case of on- 
resonance driving {5j = 0) and identical Rabi frequencies 
flj = i2. It is convenient to work in the interaction picture, 
defined by 



J = 1.2 / V J = l,2 

(8) 

Under this transformation, crj is left invariant while cr^ 
transforms as 

exp ^ ^/^crji^ crj exp ^—-Hajt j — ay cos Ht—aj sin i7t, 

and we obtain a master equation with time-dependent dis- 
sipation. 



J = l,2 

-cTjPcrJ + ct/p + pcr+ - a~p - paj) 

+ ^ E ^1 Sin(2l2t) (cr|pf7- + cr-pcr| - CT|pcr+ - CF+ pa 

J = l,2 



+ — ^ F.j[?, + Acos{m)+cos{2m)\ 

X {2ajpa+ - a+a-p - pa+aj) 

- ^ r,[3-4cos(/2t) +cos(2/?t)] 
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i ^ r,[l-cos{2nt)] {a+pa] 
ajpa- + a] pa] ~ p) . 



Wc employ now the secular approximation, il 3> Fj 
[TT] , which allows us to eliminate all the oscillating terms 
in the dissipative part fi (see Sec. 3 in 1 ,), and we obtain 

p « [2ala^^ + alal + ala^.p] + 



- {a+pa++ajpa] 



a, per. 



where F{t) = \z{t)\ - .Ja{t)d{t) and G{t) = - 
^Jh{t)c{t). In the following sections, we will study the con- 
currence evolution for some specific initial X states. To 
clearly see the limits of the secular approximation, for all 
the plots below we have taken Q = 25r, which does not 
strongly satisfy f2 ^ F. Otherwise, for larger values of i7, 
the secular approximation gives results indistinguishable 
from numerical simulations. 

From now on, we introduce the standard (exponential) 
decay amplitude r]{t) = exp{—rt/2), < 77 < 1, a nota- 
tion which will be used throughout the remaining sections 
of the paper. 



3 Werner states 



per - (if) . (10) The Werner state is defined as 



As already mentioned above, we will study the so- 
called "X" -states, which have density matrices in the form 



ps(0)=p./(0) = p(0) 



uq Wo ' 
60 20 
zo* Co 
.w^ do 



(11) 



A first observation is that this form is not preserved by 
evolution with Eq. ([5]), due to the existence of driving 
fields; thus the secular approximation brings in a welcome 
simplification, since it eliminates these fields; indeed, the 
evolution under Eq. (jlO[) preserves the density matrix in 
the "X"form. 

We can then take for the structure of the density ma- 
trix at any time t, 



Pit) 



a{t) w{t)\ 
b{t) z{t) 
z*{t) c{t) 
,w*(<) d{t) j 



(12) 



and substitute it into (fTO|) to obtain the following kinetic 
equation for the non-zero density matrix elements: 

^a{t) = \r{h{t) + c(t) - 2a{t)\ -f \u:,{w{fy - z«*(i)], 



dt 

d , , 

d 



^rKi)-26(t) + d(t)] + ^c^,[z(i)-z*(i)], 
\r{a{t)^2c{t)^d{t)\-)u:,{z(t)-z*{fyl 
^-r\h(t) + c{t) - 2d(t)\ - U^A^ify - w*{t)\ 



^^z{t) = lr[w{t) + w*&t) - 10z{t)] + ^Lu,[b{t) - c(t)], 



dt 



This system of equations can be analytically solved 
for arbitrary initial X states, and, the concurrence can be 
computed by a simplified formula [4] 



C[p](i) -2 max {0,F(t),G(t)}, 



(13) 



Pw 



1-/ 



4/-1 



3 3 

/I-/ 
(l + 2/)/2(l-4/)/2 
(l-4/)/2(l + 2/)/2 

V 







1-/. 



(14) 



where \^-) = 1/V2(|01) - |10)) is the Bell singlet state. 
The parameter / is called fidelity and takes values between 
1/4 < / < 1. Werner states with fidehty / < 1/2 are 
separable; for / > 1/2 they have a concurence of 2/ — 1. 
At zero driving and zero coupling, it is known that there is 
sudden death only in the region / < 0.714 and exponential 
decay for / > 0.714 [J. In the case of finite driving, it has 
been shown that Werner states will all decay by sudden 
death, as a consequence of the fact that driving modifies 
the structure of the decay Liouvillean . As we will see 
below, also in the case in which both the coupling and the 
driving are non-zero, this property remains valid, and all 
the states decay by sudden death. 

For the Werner state, Eq. (fT5|) has the solution 

«(0 = rf(i)-^[3+(l-4/)r;3(i)], 



m = <t) 
z{t) 
wit) 



12 

4/ 



12 

^[3 + (4/-l)r;3(t)], 

V(0 [1 + ^)1, 
- 1 



12 



-r^'it) [l-vit)]. 



This solution is independent of the qubit-qubit coupling 
LUc, which means that the enhancement of entanglement 
sudden death is also independent of lUc- This is confirmed 
numerically in Fig. [U where we present also a typical time 
evolution with finite qubit-qubit coupling for these states. 
We then obtain the following results for the expression 
appearing in the formula for the concurrence of X-states, 
Eq. inD: 



Fit) 
Git) 



4/-1 
6 

4/- 



vHt) 

1 



v\t) 



4/^1 
12 
1 

4' 



v'it) 



(15) 
(16) 
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Eq. (jl6p is negative for all values of /, therefore the 
concurrence is determined by Eq. (fT5|) : all the states with 
/ > 1/2 display sudden death due to the fact that the 
equation F{t) — has one real root. For example, the 
maximally entangled state / = 1 has a sudden death time 
ts « 0.84/r. 




co,/r 

Fig. 1. (color online) (a) Time evolution of the concurrence for 
the Werner state, with coupling Uc — 5r. (b) Time of sudden 
death of entanglement as a function of the Werner parameter 
/ of the initial state and the coupling ujc for J7 — 25r. This 
plot demonstrates numerically the analytical result that the 
time-evolution of Werner states does not depend on coupling. 
For both simulations, il — 25r. 



coupling all the YE states decay by sudden death T , and 
this feature is preserved also in the case of both nonzero 
coupling and driving, as we will see below. 

With this initial state, the solution of Eq. psp is 



ait) 

d{t) 

b{t) 
wit) 
with 



^\3~rjHt)~2{2a-l) 



cosh ( — 
4 



sinh 



Ct\ r 

C 



v'it) , 



12 



3-r/'^(t) + 2(2a-l) 



cosh 



sinh 



Ct\ r 

C 



-,-rt 



cW-^ [3 + ^7^(0], z{t)^\v\t)[l + ^{t)]. 



-(2«-l)^sinh(^^ 



1 



{I -flit)) 



(18) 



In the limit of ujc S> -T, we can approximate C w 2iL0c 
and we find 



G{t) 



' [3 - Ti^{t)Y - 4(2a - l)2ry2(i) cos2 



LJct 



vHt) + vHt) 



6 



'(2a- 1)2 sin^ 



UJct 



3 + r?3(t) 
12 



(19) 



(20) 



We note that G{t) is negative for all values of a, and 
therefore the equation F{t) = sets the value of tESD- 
A typical example of evolution is shown in Fig. (5] These 
states start with nonzero entanglement, which decays to 
zero in an time t^sD shown in the parameter space s — lUc 
in Fig. O However, unlike the diagonal states that we will 
analyze below, the YE state displays only very mild de- 
layed birth of of entanglement phenomena at very large 
LUc/r. For moderate values of ujc/T, the transitions re- 
sulting in the lobs of Fig. [3] are smooth. 



4 YE state 

The YE states have been introduced in [2], 



Pye = t; 



/I 



V 



-a 0~ 
110 
110 
a, 



(17) 



where the parameter a takes values in the interval [0, 1]. 

In the case of zero coupling and zero driving it has 
been shown [21 that only the states with a > 1/3 suf- 
fer sudden death on entanglement, while the states with 
a < 1/3 decay exponentially. For finite driving and zero 



5 Diagonal states 

An important class of states are the diagonal states, of 
which we consider two subclasses, namely states that are 
diagonal on a subspace of single excitations {i.e. only one 
quanta is present initially in the system) , and states which 
are diagonal in the subspace with zero and two excitations. 
Since g and e are typically used as indices for ground, re- 
spectively excited states, we will call Peg-ge the first class 
of states and Pee-gg the second class. Both these states 
have zero initial entanglement: as we will see however, a 
finite concurrence is acquired during the short-term evo- 
lution and decays in a finite time afterwards. 
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Fig. 2. (color online) Time evolution of the concurrence for 
the YE state, with coupling ujc = S-T and i7 = 25r. 
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Fig. 3. (color online) (a) Time of sudden death of entangle- 
ment of YE states as a function of the Werner parameter a 
of the initial state and the coupling ujc for SI — 25/', obtained 
numerically, (b) The corresponding analytical result of (a) in 
the secular approximation. 



5.1 Diagonal one-excitation states {eg — ge states) 

The generic form of eg — ge states is 



Pe3-9e = (l-p)|10)(10)+p|01)(01| = 



'0 0^ 
1 

p 

,0 0, 



(21) 

where the parameter p g [0, 1]. The physical picture that 
these states correspond to is the follovifing: suppose we 
have two identical qubits, each with its own source of dis- 
sipation. We excite one of them by a fast tt pulse (fast 
compared to F) and then we allow them to interact with 
a coupling strength lOc- For example, with quantum su- 
perconducting circuits this type of initial state is realized 
in most experiments [5J[n] . 

In the absence of dissipation we would expect the trans- 
ferring back and forth of a single quanta of excitation be- 
tween the two qubits, with periodic creation of entangle- 
ment due to coupling. The dynamics of the system would 
explore only the subspace spanned by |01) and |10), having 
at times components on the maximally entangled Bell sin- 
glet states \^^) = (|01) ± |10))/\/2. As we wiU see below, 
in the presence of dissipation and driving, this periodic 
process is still present, and if it is faster than the decay 
rate it leads to periodic creation and extinction of entan- 
glement, as the concurrence competes with the effect of 
the environmental noise. Numerical simulations are pre- 
sented in Fig. m 

For eg — ge states, the solution of Eq. (fT3|) is 





+ 2 


(!) 


r' 




i. 




- 2 


(!) 


r' 




1. 



cosh 



cosh — 



with ^ = r'^ — 36{uJc)^ « Gii^c, for ujc ^ F. 

Since now w = 0, the concurrence simplifies to 



C{p) = 2max{0,F(i)}, 



(22) 



where 



(2p - 1) sin 



-[l-ry3(t)], (23) 



which is shown in Fig. [5]Jb). We notice here that the ini- 
tial state has zero entanglement: later, some entanglement 
builds up, decays, and revives after some time. From Eq. 
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(|23|) we see that there is indeed a periodic process (show- 
ing up in the sine function) which, if it is fast enough, 
drives the function F to positive values. The period of this 
process is independent on p and has the value 2tt/3uJc, 
which agrees very well with the period of the maxima 
extracted from numerical simulations such as Fig. |4l Of 
course, in the absence of dissipation and driving, the con- 
currence reduces to the expected purely oscillatory value 
C ^{2p~ l)|sin(3wcV2)|- 







Fig. 4. (color online) (a) Time evolution of the concurrence 
for the eg — ge state, with coupling Wc ~ 5r and il = 25r. (b) 
Time evolution for ujc = lOF and O — 25r. 



Fig. [5] shows some unusual features, seen as sudden 
changes of color along the axis lOc, which mark the appear- 
ance of another lob structure. When calculating tssD, the 
code evolves the system in time, and records the last time 
value for which the concurrence was positive. We have 
verified numerically that the sudden shifts in tssD cor- 
respond to the appearance of a new entanglement revival 
(of the type shown in Fig. 3]). The evolution within the 
fcth lobe in diagram Fig. [5] corresponds to k delayed en- 
tanglement sudden birth events (for example, the value 
LUc = lO-T and f2 = 25r as in Fig.dfb) corresponds to the 
fourth lobe in Fig.[5|). 



Fig. 5. (color online) (a) Time of sudden death of entangle- 
ment for eg—ge states as a function of parameter a of the initial 
state and the coupling Uc for O = 25/^, obtained numerically, 
(b) The corresponding analytical result in the secular approx- 
imation. 



5.2 Diagonal ee — gg states 

We now consider states parametrized by s e [0,1] of the 
type 



Pee-gg - + (1 - s)|00)(00| = 



's 



,0 1 - s, 



(24) 

Physically, they generalize the situation in which two 
qubits, both in their excited state, are put in contact. 
Again, we can solve analytically Eq. (fT3|) . 



w{t) 
a{t) 



fc^sinh(^]e-" 



1 



1 + 7]^t) + 2(23 - 1) 



,z(i) = 0, 



sinh 



d(i) = -^ 1 + 77^(0-2(25-1) 



cosh 



cosh ( 
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' sinh 



where ( is defined in Eq. (|18p . This time, the concurrence 
is reduced for z{t) = to 



C{p) = 2niax{0,G(t)} 
In the limit of Wc 3> we find 



(25) 



Git) 



(1 - 2s) sin 



V'{t)-\[l-v'{t)], (26) 



which describes again an oscillatory process overlapping 
with decay. As in the analysis before, tssD will display 
revivals of entanglement with period 2'k/uJc, a statement 
that we have also checked numerically. In Fig. [6] we show 
the corresponding diagram oHesd in the s~uJc parameter 
space. As in the case of eg — ge states, the appearance of 
lobes in this diagram (sudden large Iesd) is due to the 
emergence of a new entanglement revival. For example, 
at oJc = ISF and i? = 25F as in Fig. [6l^b), we have two 
revivals and, with these parameters, we are at the second 
lobe in Fig. [71 



6 Discussion: coupled undrived systems and 
finite temperature effects 

At this point it is worth asking the following question: 
is it possible to separate the roles played by the environ- 
ment, by the driving, and by the coupling? We already 
know what happens when coupling is not present, as most 
of the literature on the topic of sudden death of entan- 
glement studies this situation. But then, if we have cou- 
pling, isn't it the case that this will produce anyway some 
form of entanglement - then, does the driving play in fact 
any role in the emergence of the structures for Iesd pre- 
sented above? The answer is yes. As pointed out already 
throughout this paper, the existence of the type of effects 
presented in this paper depends on all three ingredients: 
driving, coupling, and bipartite dissipation. Note that bi- 
partite reservoirs are indeed needed: in the case of a single 
reservoir, the decay is exponential and Iesd is infinite. 

Indeed, if one takes the master equation Eq. ([5]) and 
simulates the evolution with fi — the corresponding 
diagrams have no interesting features: there is, somewhat 
surprisingly, no dependence oHesd on ojxxi so none of the 
diagrams would be obtained (with the exception of that 
for Werner states, which has the special property of oJxx- 
independence even at finite driving). 

In the following we summarize the results obtained 
by numerical simulations for the no-driving situation. To 
make the discussion more general, we will introduce also 
temperature in the problem, 



n. 



(27) 




(b) 




I 



Fig. 6. (color online) (a) Time evolution of the concurrence 
for the ee — gg state, with coupling Uc — BF and O — 25r. (b) 
Time evolution of the concurrence for the ee — gg state, with 
coupling uic — ISr and J7 — 25r. 



(28) 



where nj are the bosonic thermal averages, 

fij = [expihajf/ksT) ~ 
and we simulate the evolution of the system under 

L^'^^ps]. (29) 



Ps 



'i\H, 



sys : 



Ps\ 



with fii = fi2 — 0. 

We start first with the observation that, in the general 
case in which we have driving, coupling, and decoherence 
in two reservoirs at finite temperature, the only effect of 
the temperature is to blur the diagrams obtained before, 
as intuitively expected. Also, for = and at zero tem- 
perature, the problem can be solved analytically [15] . Here 
we will show that finite temperature can have nontrivial 
effects, as demonstrated by the case of eg — ge states. 

For Werner states, as expected, there is no dependence 
on Wc, either at zero temperature or at finite tempera- 
ture. Also here, the only effect of temperature is that it 
decreases the sudden death time and tends to transform 
also the regions which were decaying exponentially (e.g. 
for Werner states, for 0.714 < / < 1 [3]) into sudden 
death. 
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Fig. 7. (color online) (a) Time of sudden death of entangle- 
ment for ee — gg states as a function of parameter s of the 
initial state and the coupling uoc for n = 2^r, obtained nu- 
merically, (b) The corresponding analytical result of (a) in the 
secular approximation. 

For YE states, wc find numerically that there is no 
change of concurrence as a function of Wc, either at zero 
temperature or at finite temperature. The only effect of 
finite temperature is that it makes the YE states that 
would otherwise decay exponentially at zero temperature 
(the states with 1/4 < a < 1/3) become sudden-death 
states as well. 

For ee — gg states tESD is zero, as the initial state has 
no entanglement: the coupling does not create entangle- 
ment by itself. Pumping is needed here. There is no inter- 
esting effect at finite temperature: Iesd remains zero. 

For eg—ge states at zero temperature there is no struc- 
ture in the Iesd diagram, again due to the fact that the 
states decay exponentially. Also at very large tempera- 
tures (n ^ 1) there is no structure: this time because the 
decay is too fast and Iesd becomes essentially zero. How- 
ever, for intermediate temperatures a structure somewhat 
similar with that obtained in the case of driving emerges, 
shown in Fig. [8] (b) for ni = n2 = 0.25. 

This is somehow unexpected: usually the situation is 
that features appearing at a certain temperature would be 
enhanced at even lower temperatures (or in other words 
the only effect of temperature is usually to wash out the 
features seen at zero or very low temperatures). We note 
here the similarity between Fig. [5] (b) and the case of zero- 
temperature driving Fig. [5] (although differences in the 



values of tssD and in the period of the lobes do exist). 
This can be qualitatively understood by using the fact 
that driving the system is mathematically equivalent with 
a combination of temperature, dephasing, and squeezing 
[7j, and therefore, for certain initial states, the driving 
could lead to similar effect as temperature. 

The reason for the appearance of the lobs is then sim- 
ilar to that given in Section 15.11 consider for example 
p — 0, which corresponds to the case in which one qubit 
in the excited state is connected to the other qubit ly- 
ing in the ground state. What happens is that the sin- 
gle quanta of energy start to oscillate between the two 
qubits. The dipole-dipole form of the interaction assures 
that the dynamics occurs only in the subspace spanned 
by |01) and |10), and, during its evolution, the state of 
the system goes through Bell entangled states of the type 
|iZ/±) (|01) ± |10))/V2. Now, if there is no dissipation, 
the concurrence will simply oscillate. If the temperature is 
zero, the concurrence decays exponentially. Once tempera- 
ture is introduced, the system will display sudden death of 
entanglement and will acquire a finite Iesd- Moreover, de- 
layed sudden birth of entanglement is produced due to the 
fact that during oscillations the concurrence could go be- 
low zero but then revive after some time, winning against 
the thermal noise. 



7 Conclusion 

We have investigated the phenomenon of sudden death of 
entanglement for certain subclasses of so-called X states 
under continuous driving and coupling. We show that each 
of these states present notable peculiarities when analyz- 
ing the time evolution: for Werner states, the evolution is 
independent of coupling, for YE states the time of sud- 
den death presents soft oscillations as a result of coupling, 
and for diagonal {eg — ge and ee — gg) states Iesd shows 
abrupt oscillations, due to the sudden appearance of de- 
layed birth of entanglement for certain value of the cou- 
pling. We present both analytical formulae and numerical 
evidence for these phenomena. 
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